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Simulation of the magnetic relaxation in a model of hard superconductor reveals a new time scale below the 
irreversibility temperature. The relaxation in this new time scale, which appears in an intermediate time window, 
is a power law and is related to self organization of the magnetic flux around the critical current density J c during 
the relaxation. Emergence of a new time scale which along with the time scale for long time relaxation due to 
thermally activated process dynamically identifies the irreversibility temperature. 



1. Introduction 

Magnetic relaxation in high-T c supercon- 
ductors is experimentally known to be non- 
logarithmic over a large time interval jl| . Though 
the non-logarithmic behaviour is related to a di- 
verging pinning potential within vortex glass the- 
ory H and collective creep theory ||, the details 
of the experimental behaviour is not consistent 
with the theoretical predictions. Specifically, the 
normalized relaxation rate is found to be tem- 
perature independent with a value that cannot 
be accounted by the theory Moreover, no 
clear picture exists for the change in relaxation 
behaviour across the irreversibility temperature 
Ti rr and into a true superconducting state with 
finite persistent current. 



2. The model 

In order to understand generic processes in- 
volved during the relaxation in hard type-II 
superconductors, we have simulated relaxation 
of the thermoremanent magnetisation in a 2D 
Josephson junction array at finite temperatures. 
As have been shown previously Q] , a 2D Joseph- 
son junction array (JJA) with finite screening cur- 
rent is a realistic model for magnetic behaviour of 
hard superconductors. We consider a 2D JJA of 
size N X N in a transverse magnetic field. The 
equation of motion for the gauge invariant phase 



difference across a junction (f> is given by 
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Here, the current in a cell /,„ is scaled by the crit- 
ical current I c of the junction. The r 

represents dimensionless time, and X 2 } 
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is the dimensionless penetration depth analogous 
to London's penetration depth of a bulk super- 
conductor (<l>o is a quantum of flux, Lo is the 
self-inductance of a unit cell, and R is the nor- 
mal state resistance of the junction). The applied 
magnetic flux is represented by / = & ext /&o- The 
matrix M is directed loop-sum operator and is 
equivalent to lattice curl operation. The temper- 
ature bath is attached through the noise term 
X{t) with (X r {t)) = and (X r (r)X r .(r')> = 
2T8(t — r')(5r ir '. The temperature is in units 
of / c $o/27rfcs. We employed free-end boundary 
conditions for all the variables. Further details 
regarding the simulation and setting up of the 
equation can be found in ref . || . 

3. Results and discussions 

We present here the results for an array of 
size N — 16. The Fig.?? shows the relax- 
ation over 6 decades in r for some selected tem- 
peratures. The relaxation behaviour changes 
markedly across two temperatures : T cr « 0.24 



and T sc w 0.04. At T cr , the relaxation curve 
develops a kink in an intermediate time inter- 
val which at lower temperature develops into a 
plateau on which the dynamics is almost frozen. 
The width of this plateau thus sets a new time 
scale Tp for the relaxation of the magnetisation. 
Experimentally, the Ti rr is obtained as the tem- 
perature at which the field cooled (FC) and zero 
field cooled (ZFC) susceptibility xCO differs. In- 
set of Fig.?? shows x(T) for the model. The tem- 
perature T cr at which the kink appears in the 
relaxation curve is also the temperature at which 
the | Xfc ~Xzfc |> 0. This allow us to conclude 
that T cr is the irreversibility temperature appear- 
ing in the relaxation. This is also consistent with 
the previous simulation study in which magnetic 
irreversibility is found to set below this temper- 
ature ||. For T sc < T < T cr , the magnetisation 
M(t — > oo ) = 0. At T sc , the relaxation is frozen 
on the plateau as r — > oo, thus establishing a fi- 
nite persistent current density (hence, remanent 
magnetisation). We consider this as the transi- 
tion into the true superconducting state. By ex- 
plicitly observing the flux distribution, we find 
that the plateau is related to the crossover from 
a supercritical state with J > J c to a subcritical 
state with J < J c . We attribute this crossover 
as arising due to self-organization during the re- 
laxation process. The Tp increases rapidly with 
decreasing T. 

For T sc < T, the long time relaxation fits to 
exp[— (t/t q ) q ]. The characteristic time scale r a 
increases rapidly below T cr . One also observes 
log(r/ro) behaviour for 1-2 decades far from the 
plateau. This is the regime in which thermal ac- 
tivation is a dominant process and can be used 
to extract the characteristic time scale To for re- 
laxation at long time. Fig.?? inset (b) shows the 
temperature dependence of t^ 1 oc (magnetic dif- 
fusivity) and is plotted as a function of T — T sc . 
Though To follows the Arrhenius law, a power law 
at low temperatures as indicated by the straight 
line in the plot is surprising. The T sc obtained 
from fitting the time scale for long time relax- 
ation matches well with the T at which M(r — > 
oo) = Mq > 0. The T cr appears as the temper- 
ature at which Tq -1 deviates from the power law 
as marked in Fig. 1 inset(b). 



In conclusion, we have shown that a new time 
scale govern the flux dynamics below Ti rr which 
arises due to self-organization of the magnetic flux 
during the relaxation. This conclusively proves 
that Ti rr is associated with the crossover in dy- 
namical behaviour for magnetic flux in type-II su- 
perconductors. 

Figure Caption The normalized magnetisa- 
tion M(t)/M(0) for some selected temperatures 
(marked along the curves) for / = 5 and N = 16. 
Inset: (a)shows the x(T) for FC and ZFC condi- 
tions (T^r is marked on the curve), (b) The t^ 1 
and Tg as a function of T—T sc on log-log plot. The 
vertical axis label indicates the exponent only. 
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